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C^i Abstract 

-i— > ■ 

We analyze the roles of initial correlations between the two-qubit system and a dissipative cavity on 



quantum discord dynamics of two qubits. Considering two initial system-cavity states, we show that the 
initial system-cavity correlations not only can initially increase the two-qubit quantum discord but also 



would lead to a larger long-time quantum discord asymptotic value. Moreover, quantum discord due to 

(N ! 

initial correlations is more robust than the case of the initial factorized state. Finally, we show the initial 
^" ■ 

correlations' importance for dynamics behaviors of mutual information and classical correlation. 

o. 

O . PACS numbers: 03.65.Ta, 03.65.Yz, 03.67.Mn, 05.70.Fh 



^ Introduction: Recently, lots of interest have introduced what they called quantum discord, 
been devoted to the definition and understand- They marked the beginning of a new line of re- 
ing of correlations in the quantum systems. En- search shifting the attention from the entangle- 
tanglement is a kind of quantum correlation that ment vs. separability dichotomy to the quantum 
has been playing a central role in quantum infor- vs. classical paradigm. The quantum discord 
mation and communication theory [1]. However measures quantum correlations of a more gen- 
there are other nonclassical correlations apart eral type than entanglement, there exists sepa- 
from entanglement [2-4] that can be of great im- rable mixed states having nonzero discord [5]. 
portance to these fields. In order to characterize Interestingly, it has been proven both theoret- 
all nonclassical correlations, Ollivier and Zurek ically and experimentally that such states pro- 

vide computational speedup compared to classi- 
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cal states in some quantum computation models necessary. 

[5,6]. The influence of initial correlations on the 
Understanding of the quantum dynamics of open system dynamics has recently been inten- 
open systems is a very important task in many sively studied [16-22]. In order to clear what 
areas of physics ranging from quantum optics happens to the quantum discord in this situa- 
to quantum information processing and to quan- tion, we study an exactly solvable model for the 
turn cosmology. Therefore, the study of quan- time evolution of two atoms interacting with a 
turn and classical correlations' dynamical be- lossy cavity. In this model where the coupling 
haviors in the presence of both Markovian [7,8] strength effects are considerable, the presence 
and non-Markovian [9,10] decoherence has at- of system-cavity correlations invalidates the ini- 
tracted much attention in recent years. It is be- tial state in which system and cavity are inde- 
lieved that the quantum correlations measured pendent. We investigate quantum discord dy- 
by the quantum discord, in the Markovian case, namics for two initial states of two atoms and 
decay exponentially in time and vanish only the lossy cavity, in which the initial reduced 
asymptotically [11], contrary to the entangle- density matrixes of the atomic system and the 
ment dynamics where sudden death may occur cavity are the same. Under a certain condition, 
[12]. In particular, Refs. [13,14] have discov- the initial system-cavity correlations not only 
ered that quantum discord can be completely can initially increase the atomic quantum dis- 
unaffected by Markovian depolarizing chan- cord but also would lead to a larger long-time 
nels or non-Markovian depolarizing channels quantum discord asymptotic value. For another 
for long intervals of time, and this phenomenon condition, correlations between two atoms and 
has been observed experimentally [15]. In these the lossy cavity, can more effectively restrain 
above studies, quantum discord dynamics of the reduction of quantum discord than the case 
open system in which system and environment of the factorized state. And then we also ana- 
are initially separable has been analyzed in de- lyze the different dynamics behaviors of mutual 
tail. As everyone knows it is possible, and of- information and classical correlation due to the 
ten unavoidable in experiment, to create corre- initial system-cavity correlations. These find- 
lations between system and environment, there- ings obtained from our article show that, dy- 
fore such correlations will play important roles namics of mutual information, quantum discord 
in the time evolution of the system. So the study and classical correlation not only depend on the 
of quantum correlation dynamics due to the ini- system degrees of freedom, the initial system- 
tial system-environment correlations is certainly cavity correlations must also be properly taken 



into account. we measure the subsystem Y by a complete set 
Quantum discord and classical correlation: We of projectors JT i? corresponding to the outcome 
now present a brief review of the classical i, which yields p x \i = Tr Y (Y[iPxY Yli)/Pu 
correlation and quantum discord. In classi- with p ri = Tr XY (Y[iPxYY[i)- Then the quan- 
cal information theory, the information can turn mutual information can alternatively de- 
be quantified by Shannon entropy H(X) = fined by 

-J2 x P\x= x logP lx=x , where P\ x=x is the prob- ^ {x . y) = 5( ^ } _ g^y^ (3) 

ability with X being x. Similarly, the joint en- 



tropy, which measures the total uncertainty of 
a pair of random variables X and Y, is defined 

as H(X, Y) = - J2 x , y P\X= X ,Y=yl0gP\ X = X ,Y=y, 

with P\ x =x,Y=y being the probability in the case 
of X = x and Y — y. Then the total cor- 
relation between X and Y can be measured 



where S U .(X\Y) = J2iPi S (px\i) is condi- 
tional entropy of the quantum state. The above 
quantity strongly depends on the choice of the 
measurements {n?}- By maximizing Jy\.(X : 
Y) over all Uli}' we define the classical corre- 
lation between X and Y 



by the mutual information which is defined as C(X : Y) = max Jtt. (X : Y), (4) 

I(X : Y) = H(X) + H(Y) - H(X, Y), whose and ^ qumtum ^ as 

quantum version can be written as [23] 

V(X :Y)=1(X :Y)-C(X :Y), (5) 

l(X:Y) = S(p x ) + S(p Y )-S(p XY ), (1) 

which is interpreted as a measure of the quan- 

where S(p) = —Tr(p\ogp) is the von Neu- . , .. ro A1 T( ^ . , f 

v ' ' y ' t> ' ' turn correlation |2-4J. It is zero only for states 

mann entropy of p, and Px(py) is the reduced ... , . , , .. , f 

^ J ' r vr ' with classical correlations and nonzero for states 

density matrix of p X y by tracing out Y(X). By V1 , . T ^ , 

j rAi j & v j j wl th quantum correlations. In particular, quan- 

introducing the conditional entropy H(X \Y) = ,■ ,■ , . ., , , , rr 

° trj v ' ' turn discord is equal to the entanglement of for- 

H(X,Y) — H(Y), we can rewrite the mutual <, **•..■ , . r . 

v ' ' v ' mation for pure states, it is not true for mixed 

information as . . . . . a -. 

states, since some states present finite quantum 

I(X : Y) = H(X) - H(X\Y), (2) discord even without entanglement [2]. 

, rr/v , v v ^ TTlv ,^ r s The model: We consider two atoms A and B in- 
where H(X\Y) = Y. y PY= y H(X\Y = y) = 

v^ , . „, ,.„. , teracting with a dissipative cavity. The Hamil- 
-l^x,yPx=x,Y=ylospx=x,Y=y is the conditional 

,. , , . , , v avc ax. tonian of such total system in the rotating-wave 

entropy of the random variables X and Y for the J ° 



approximation is given by H = H + H,, 



int i 



average uncertainty about the value of X given 

., ., . fV . , , which, in the basis { I ##), leg), lye), lee)} reads 

that the value of Y is known. In order to gener- u ^ J/ ' J/ ' •* " ' /J 

alize the above equation to the quantum domain, H = u (a A a A + cr+cr B ) + u c a^a, (6) 



H int = tt(a+a + o+a) + h.c, (7) tern contains at most two excitations. In this 

case the dynamics of two atoms can be effec- 
tively described by a four-state system in which 



. 4 b ^ i ja t. t- ■ ■ case the dynamics of two atoms can be effec- 

here a±, a± are, respectively, the Pauli raising J 



and lowering operators for atoms A and B, u>q is 

^ r, , r r <_ j + three states are coupled to the cavity mode in a 

the Bohr frequency of two atoms, a and a ] are r J 

^. .. ., , r ., ladder configuration, and one state is completely 

the annihilation and creation operators tor the ° r J 

, , . , . , . , , ^ j, decoupled from the other states and form the 

cavity mode, which is characterized by the fre- r 

, . ,. „ „ . field. In the basis {10) = \gg),\+) = (leg) + 

quency u c and the coupling constant \l. For the u/ i-^/'i i m ^/ 

v f • r v • t u f n • lffe))/\/2,|-) = (|ey)-|5e))/V^,|2) = |ee)}, 

sake of simplicity, in the following we assume ,a '" ' ' Vl ' '^ /y/ ' ' ' /J 

., . ., ., V1 ., the Hamiltonians (6) and (7) can be rewritten 
that the two atoms interact resonantly with the 

dissipative cavity mode. ff , = 2wo | 2)(2 | +a;o( | +)(+ | + |_ )( _| )+ ^ a;fea t afc5 

In this section, we investigate the dynamics fc 

(9) 

of the two atoms interacting with a dissipative ^ 

H' int = V2tt(\+)(0\a+\2)(+\a) + h.c. (10) 
cavity by making use of the master equation. 

Firstly, we focus on the case in which the total From tne total Hamiltonian given by © and 

system contains only one excitation. In the in- CUB, the subradiant state |-) does not decay, 

teraction picture, the exact dynamics of the two and tne superradiant state is coupled to states 

atoms is contained in the following master equa- 1 0) and 1 2) via the cavity mode. The transitions 

(i on |0) — > |+) and |+) — > 1 2) have the same fre- 

^ Y quencies and identically coupled with the cav- 

— = -i[H int ,p] - -[a ] ap - 2apa ] + pa) <i$) , . 

at 2 lty. In the interaction picture, the dynamics of 

where p is the density operator for the two atoms two atoms interacting with a lossy cavity can be 

and the cavity mode, V is the decay rate of the treated in the following master equation 

cavity mode. In order to find the dynamics of dp Y 

-17 = - * [Hint, P] - 77 W a P ~ 2apa T + Mv-) 

two atoms, we solve the master equation in Eq. at A 

(8). In the case which the total system con- Focusing on the case in which there are at most 

tains only one excitation, we need to solve a two excitations in the total system, we should 

set of 16 differential equations obtained from solve a set of 64 differential equations obtained 

Eq. (8). Then, tracing out the cavity degree of from Eq. (12). The symmetry properties of the 

freedom, we obtain the reduced density matrix system allow to transform the set of differential 

of the atomic system for the total system which equations into decoupled subset of differential 

only contains one excitation. equations of smaller size. Finally, we obtain the 

Then we treat another case that the total sys- density matrix of the reduced atomic system by 



tracing out the cavity degree of freedom. the total system contains only one excitation. 

To illustrate the roles of the initial correla- So the density matrix of the total system at time 
tions between two atoms and cavity on the quan- t can be calculated by Eq. (8). However, not 
turn discord dynamics of two atoms, we con- only the basis (|<?e0), |e<?0) and \ggl)) but also 
sider a initial condition p ABc = \^)ab c {^\, with those basis (\gel) and \egl)) are contained in 
\^)abc = o:\geO) + (3\eg0) + j\ggl) (having the factorized state, so p ABc actually consists 
correlations), and here \a\ 2 + \(3\ 2 + \^\ 2 = 1 of two parts, one satisfying one excitation, 
(a,/3,7 7^ 0), |e) and \g) are the exited state and the other meeting two excitations. Hence 
and ground state of atoms, |0) and |1) are the the evolutional density matrix of p ABc can be 
vacuum state and the single-photon state of the acquired by Eqs. (8) and (12). Tracing out 
lossy cavity. Obviously, two atoms and the cav- the cavity degree of freedom, we obtain the 
ity in this initial state is density matrix of the reduced atomic system 

for these two different initial states. In the 
basis {\gg), \eg), \ge), |ee)}, we measure the 
(12) atom B from the matrix pAB(t) by projecting 
on {cos-#|e) B + e^sim!%) B , e -z ^sini?|e).B — 
cos i?|g , )s}. Then the quantum discord and 



Pab 



( \i\ 2 o\ 

\a\ 2 a(3* 

a*/3 \/3\ 2 

\ 0/ 



Ul 2 _|_ I «|2 q \ classical correlation could be calculated numer- 

Pc = I • (13) 

Ul 2 / ically using Eqs. (4) and (5). 

Therefore, the initial correlations must illustrate Numerical results and discussions: In this 

clear effects on the dynamics of the atomic article we stud y a s y stem whose d y nam " 

quantum discord. In the following, by com- ics is described by the well-known damped 

paring to the second initial condition that the Tavis-Cummings model. To see the effects 

two-atom system and the cavity are in the of initial system-cavity correlations explicitly, 

factorized state p% c = p AB <g> p c , we mainly We consider two distinct initial states, the 

study the different atomic quantum discord system-cavity correlated state p« c and the 

dynamical behaviors due to these two initial factorized state p% c . The reduced density 

atoms-cavity states, which contain the identical matrices for both the atomic system and the 

states of subsystems. loss y cavit y ** the same - For simplicity, 

For p% c , it is composed by the basis |<?e0), we choose the parameters a = sin 9 cos <p, 

\eg0) and \ggl), that is to say the initial corre- f 3 = sin ° sin V^ and 7 = cos 6 (6 e [0, tt/2] and 

lation state Pa % c satisfies the first case in which f G t ' *"])• In Fl 2- l > b y choosing p {1 A \ c and 
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FIG. 1 . Time evolution of quantum discord of two 
atoms as a function of the dimensionless quantity Qt 
with 9 = 7r/3 and ip = 3tt/A (a) the initial system- 
cavity correlated state p AB as the total initial state, 

(2) 

(b) the factorized state p ABc as the total initial state. 

(2) 

Pabc as ^ initial states of the total system, we 
plot the time evolution of quantum discord for 
two atoms as a function of the dimensionless 
quantity Vtt, with 9 = n/3 and ip = 37r/4. We 
can observe that the atomic quantum discord 
can be initially increased to a maximum in 
the case p ABc (as shown in Fig. 1(a)), which 
contains system-cavity correlations at t — 0. 
While in the factorized initial state p ABc , the 
initially decreased atomic quantum discord to a 
minimum would happen (as shown in Fig. 1(b)). 
Hereafter, how much the coupling strength of 
our regime influences the quantum discord? 



Since T is connected to the lossy cavity corre- 
lation time r B by the relation r B w 1/T, and 
il is related to the time scale t r over which 
the state of the system changes, r R w I/O,. 
So T/Q < 2 characterizes the strong coupling 
regime, and T/Cl > 2 means the weak coupling 
regime. Fig. la presents the case for the initial 
system-cavity correlated state p ABc , it can be 
seen that the initially increased quantum discord 
maximum increases with the decrease of the 
T/fl corresponding to the enhancement of the 
coupling strength. We also clearly see that for 
the factorized initial state p ABc in Fig. lb, the 
initially decreased quantum discord minimum 
decreases with the decrease of the F/Q. 

Then, by investigating the subsystem of two 
atoms, we shall compare the roles of different 
initial system-cavity states (p ABc and p ABc ) 
on the atomic quantum discord dynamics. 
Although the initial reduced density matrices 
for both the atomic system and the lossy cavity 
are the same as in the full calculation, Fig. 2 
shows that the presence of the system-cavity 
correlations in the initial state changes the 
atomic quantum discord dynamics dramati- 
cally. In figs. 2(a) and 2(b), quantum discord 
dynamics is given for the conditions 9 = n/3 
and if = 37r/4 in the strong coupling regime 
and weak coupling regime, respectively. For 
the case p ABc , the atomic quantum discord can 
firstly increase to a maximum, then periodically 
decrease to a long-time asymptotic value in 



the strong coupling regime. While in the weak 
coupling regime, quantum discord would reach 
a maximum at first and then gradually decrease 
to the asymptotic fixed value. For the other case 
Pabc quantum discord of two atoms initially re- 
duce to a minimum, and eventually periodically 
increase to another long-time asymptotic value 
in the strong coupling regime (no oscillations 
are present in the weak coupling regime). Due 
to the subradiant state |— ) in two initial states 
does not decay in the evolution process, the 
long-time asymptotic quantum discord would 
be acquired. The final long-time asymptotic 
value due to the system-cavity correlations is 
much larger than the case of the factorized 
state. So we can conclude that the initial 
system-cavity correlations not only can initially 
increase the atomic quantum discord but also 
would lead to a larger long-time quantum 
discord asymptotic value. 

Figs. 2(c) and 2(d) give the atomic quantum 
discord dynamics for the conditions 9 = n/3 
and ip = 7r/4 in the strong coupling regime 
and weak coupling regime, respectively. 
When if = tt/4, the initial atomic state 
does not contain the subradiant state |— ), the 
asymptotic stationary state of two atoms is 
pAB(t -> oo) = \gg){gg\, so the long-time 
asymptotic value of the atomic quantum discord 
is zero. For both the initial system-cavity 
correlated state or the factorized state, in the 
strong coupling regime the atomic quantum 



strong coupling regime (r=0.2ii) 




.'•n«/»ii.'-- f " v " 
lui-'- 



weak coupling regime U" 5£l) 




FIG. 2. Time evolution of quantum discord of two 
atoms as a function of the dimensionless quantity Qt 
with 6 = it/3, (a) if = 3vr/4 and T = 0.2ft, (b) 
ip = 3vr/4 and F = 5ft, (c) ip = vr/4 and T = 0.2ft 
and (d) ip = it /A and Y = 5ft. For the cases of 
(i) the initial system-cavity correlated state p ABc ( 

(2) 

red solid curve), (ii) the factorized state p ABc (dark 
dash-dotted curve). 

discord presents damped oscillations while in 
the weak coupling regime quantum discord 
decays only asymptotically to zero. A compar- 
ison between the red solid curve and the dark 
dash-dotted curve in Figs. 2(c) and 2(d) reveals 
that for the same initial atomic state, the atomic 
quantum discord dynamics due to the initial 
system-cavity correlations is more robust than 
the case of the factorized state. This suggests 
that correlations between two atoms and the 
lossy cavity, effectively restrains the reduction 
of atomic quantum discord. 
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FIG. 3. Correlation (mutual information, quan- 
tum discord, classical correlation) dynamics of two 
atoms as a function of the dimensionless quantity 
ttt in the weak coupling regime with 9 = ir/3, (a) 
<f = 3vr/4 and T = 50, (b) 92 = vr/4 and T = 50. 
For the cases of (i) the initial system-cavity corre- 
lated state p ABc ( the solid curves a, c and e), (ii) 
the factorized state p ABc ( the dash-dotted curves b, 
d and /). 

To further understand the roles of different 
initial system-cavity states (p ABc and p ABc ) on 
dynamics of the total correlations, we study 
the atomic mutual information X and classical 
correlation C, respectively. Fig. 3 shows the time 
evolution of the atomic mutual information, 
quantum discord and classical correlation in the 
weak coupling regime. It is clearly found that 



for the initial correlated state p ABc , the atomic 
classical correlation is greater than quantum 
discord except in the initial period of time. 
While the classical correlation is always less 
than quantum discord for the initial factorized 
state p ABc . When tp = 3n/A, because of the ini- 
tial system-cavity correlations, both the atomic 
classical correlation and mutual information 
can firstly increase to a maximum, and then 
get to a fixed long-time asymptotic value. In 
addition, the fixed long-time asymptotic values 
of classical correlation and mutual information 
due to p ABc are much larger than the case of 
Pabc (Fig- 3(a)). This result is same as quantum 
discord dynamics. When cp = ir/4, the atomic 
classical correlation and mutual information 
due to the correlated state p ABc would increase 
to a maximum value at first, then reduce 
asymptotically to zero. This finding is different 
from quantum discord dynamics, which only 

decays to zero (Fig. 3(b)). Comparing to the 

(2) 
case of the initial factorized state p ABc , classical 

correlation and mutual information of two 
atoms also can be more robust. 
Conclusion: In conclusion, we have studied 
quantum discord of two qubits in the presence 
of initial system-cavity correlations. The cor- 
relations between the qubits and a lossy cavity 
is found to have important effect on the time 
evolution of quantum discord when two qubits 
couple with a common lossy cavity. Finally, we 
also analyze the different dynamics behaviors 
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of mutual information and classical correlation 
due to the initial system-cavity correlations. 
In comparison with some recent work on the 
initial factorized state between the system 
and environment, our present work might be 
more practical to explain the total correlation 
dynamics behaviors of the system. 
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